Performance bound for quantum absorption refrigerators 
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An implementation of quantum absorption ciiillers witii three qubits iias been recently proposed, that is ide- 
ally able to reach the Carnot performance regime. Here we study the working efficiency of such self-contained 
refrigerators, adopting a consistent treatment of dissipation effects. We demonstrate that the coefficient of per- 
formance at maximum cooling power is upper bounded by 3/4 of the Camot performance. The result is inde- 
pendent of the details of the system and the equilibrium temperatures of the external baths. We provide design 
prescriptions that saturate the bound in the limit of a large difference between the operating temperatures. Our 
study suggests that delocalized dissipation, which must be taken into account for a proper modelling of the 
machine-baths interaction, is a fundamental source of irreversibility which prevents the refrigerator from ap- 
proaching the Camot performance arbitrarily closely in practice. The potential role of quantum correlations in 
the operation of these machines is also investigated. 



PACS numbers: 03.65.-w, 03.65.Yz, 05.70.Ln, 03.67.-a 

Introduction. — The study of quantum thermal machines has 
attracted an increasing attention over the last years. This is 
motivated on one hand by the fundamental interest in un- 
derstanding the emergence of basic thermodynamical princi- 
ples at the quantum mechanical level lfT HT2l . and on the other 
hand, by the potential technological applications of these ma- 
chines, for instance to control the heat transport in nanoengi- 
neered devices lfT3l - [T6l . In particular, several models have 
been proposed fpJTl realizing quantum absorption chillers, 
that is, refrigerators in which the external source of work is 
replaced by a heat bath. A very simple realization of such 
refrigerators, which has been introduced in |5 -7|, consists of 
three interacting qubits, with a vanishingly small interaction 
strength, each one in contact with a heat bath. In spite of the 
technological challenges behind its physical implementation, 
this machine can be experimentally realized, e.g., with super- 
conducting qubits or arrays of quantum dots lfT4l [T5]| . Fur- 
thermore, its operation may be understood in a very neat way, 
providing a physical insight into the sources of irreversibil- 
ity in absorption chillers Q. It has been predicted that, with 
a suitable choice of the machine parameters, such a refrig- 
erator can ideally attain a coefficient of performance (COP) 
reaching the Carnot bound eq |6|. However, we argue that the 
central assumption of vanishing mutual interaction between 
the refrigerator qubits cannot be realistically maintained. As 
long as the interaction is finite, each bath will exchange energy 
with the whole three-qubit system, rather than just locally with 
the single qubit to which it is connected. This is usually the 
case with any interacting multiparticle dissipative system. As 
we shall show, the resulting delocalized dissipation prevents 
the refrigerator from approaching the Carnot limit arbitrarily 
closely, thus embodying a fundamental source of irreversibil- 
ity that will arise in all concrete implementations. 

This situation is reminiscent to that of realistic heat engines 
or cooling cycles, topical areas of study in finite -time thermo- 



dynamics. There, the finite heat transfer rates constitute the es- 
sential source of irreversibility which makes the Carnot bound 
unattainable in practice. For this reason, an important line 
of research deals with devising alternative, tight performance 
bounds, such that some suitable figure of merit of the machine 
under consideration is maximized f T7142TII . In this spirit, we 
address the following question: What is the highest achievable 
COP at maximum cooling power for the quantum absorption 
refrigerators of Refs. fS'-T]? Answering this question would 
provide a practical performance bound against which the ef- 
ficiency of any future realization of these machines could be 
benchmarked. Another relevant and related question to ask is 
whether the "quantumness" of the refrigerator, as revealed for 
instance by the stationary quantum discord Il22ll23l . plays any 
role in its operation. This would help to unveil connections 
between quantum correlations and efficient energy transport 
out of equilibrium, that so far have remained elusive. 

In this Letter we answer both questions. In the first place, 
by considering unstructured bosonic baths and a consistent 
dissipative qubit-bath interaction, we find that the COP at 
maximum power is tightly upper bounded by 3ec/4, where 
£c = (l - 2^)/(Z^ _ is the Carnot COP and {T„, T;,, r„,) 
are the three temperatures between which the refrigerator op- 
erates. We also give sufficient conditions to saturate this 
bound in the limit of large temperature difference Tc/Ti, «: 1. 
Secondly, we issue a comprehensive analysis of stationary 
quantum correlations in the various relevant qubit-qubit parti- 
tions. Although a nonvanishing discord is always found in a 
specific partition, it does not relate with the stationary heat 
flows, reinforcing the idea that this family of thermal ma- 
chines operates in an effectively classical way [T^, despite hav- 
ing a genuinely quantum physical support. 

Preliminaries. — Let us begin by introducing the total Hamil- 
tonian for the refrigerator The Hilbert space of the system is 
•Ks = "J^s.w ® "Ks./i ® 'Ks.ci where we label the three qubits as 
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FIG. 1 : (Color online). Schematic representation of the three-qubit 
absorption refrigerator. The refrigerator qubits dissipate into their 
respective baths, with equilibrium temperatures r„ > T;, > T^, at a 
rate y. The three-body interaction, of strength g, allows for energy 
exchange between the refrigerator qubits, whose energy spacings aia 
are required to satisfy a>i, = 0)^ + aj,^.. 



'work', 'hot', and 'cold' (w, h, c), after the heat baths to which 
each of them is connected (see Fig. [T}. Their free Hamilto- 
nians are Hda - |1q.) where a - {w, h, c], and we 
work in natural units ti - kg - 1 . The corresponding bosonic 
equilibrium reservoirs are given by Hb.c - YjA ^a^I ifla-.i- As 
local qubit-bath dissipative coupling, we choose terms of the 
form HD,a = Vt('^.i„o-x„ + Cv„cryJ (gi 2i g^iaaj - al j, where 
{c^^,, CyJ e R, and oc yfojA to ensure flat spectral densi- 
ties J(a)) ~ g1^JoJa,A |24|. Here, we absorbed the order of 
magnitude of J (oj) into the dissipation rate j. With no loss 
of generality, we can set c^^ = 1 and Cv„ - 0. This kind of 
system-environment coupling stands e.g. for the dipole in- 
teraction between a two-level atom and the electromagnetic 
field at thermal equilibrium f2A\. It only remains to spec- 
ify the three body interaction between the qubits, given by 
Hi - §(|lvvO/,lc)(0,vl/,Oc| + h.c), where g is the interaction 
strength. The qubit energies are chosen as Uh - u)c + Wh 
(<jj/, > (On), so that the transition |l,vO;,lc) |0h1/i0c) is ap- 
proximately resonant, as long as g <sc 1 . The total Hamilto- 
nian is finaUy Hj = Ho,a + H, + 2„ Ho,a + 2» 

We now briefly explain the operation of the refrigerator (see 
Q for details). In the ideal scenario of vanishing g, the re- 
duced stationary state of the work and hot qubits p"^^ = Tr^ g, 
has high fidelity with t„, ® t/,, where = z^^e^^""^^" 
stands for the thermal state of qubit a at the equilibrium 
temperature r„, and g denotes the reduced state of the three 
qubits after tracing out the heat baths. It follows that the 
truncation of ^p^^^ into the two-dimensional subspace TYs ,, 
of 'Hs^w ® 'J^sj: spanned by {11^,0/,) , |0„1/,)), which defines a 
'machine virtual qubit' y, has an effective virtual temperature 
approximately given by T,, = (ojh - w,,) / (tOh/Th - io^./T^). 
The interaction Hj allows the cold qubit to exchange energy 
with the machine virtual qubit, while being simultaneously 
thermalized by the cold bath through Hj) c- Suitable choice 
of frequencies and temperatures may result in T,, < Tc, so 
that the (non-equilibrium) stationary state g"^ is effectively 



colder than Tc- The excited state population deficit in g^ is 
compensated by a net energy transfer from the cold bath (that 
stands for the object to cool) into the refrigerator This is 
what we shall understand by cooling. The machine therefore 
just mediates between the cold object at temperature Tc and 
a suitably filtered virtual temperature Ty. Thermalization is 
then responsible for cooling within the window < T,. < Tc- 
When Ty - Tc and always under the assumption of localized 
dissipation, which is consistently realized only for vanishing 
g, the refrigerator saturates the Carnot bound sc on COP 161. 

Realistic modelling of the dissipation. — We shall now ex- 
tend the model of Refs. |5 -7| to take into account a consis- 
tent treatment of dissipation, necessary for the realistic case 
of even very small but nonvanishing g. We can derive a 
general equation of motion for the qubits (in the form g - 
-i[HKi,g] + Tja^a [g]), from first piinciples, by employing 
the standard Born-Markov assumption of weak, memoryless 
system-environment coupling 1 24| . Here i/ref = Ziu- H^) a+Hj. 
Such master equation is written as ||25l 

g = -/ [//ref, g]+^ ^a,aj {Aa.ajgAl^cj ' sMl.^^^ff.o;, £)) + ) ■ (1) 

The spectral correlation tensor, denoted by Fq.^, is propor- 
tional to the real part of the power spectra of the bath cor- 
relation functions. Note that since the heat baths are inde- 
pendent, the correlation tensor is diagonal in a. We use its 
explicit form for electromagnetic radiation at thermal equi- 
librium: r„_^ oc exp(w/27'o,) (sinh w/2rQ,)"' |24|. The 
non-Hermitian Lindblad operator Aa^^, associated with bath 
a, performs transitions of frequency o) at rate Fa^ ^j, between 
the (^-independent) eigenstates of the refrigerator Hamil- 
tonian i/ref- They result from the decomposition of the 
system-environment couplings ( yfycTx^ - 2„ Aq,,^) as eigen- 
operators of the refrigerator Hamiltonian ([//let, Aa,ai\ - 
-djAa^iS)- Note that the corrections to H^^f resulting from the 
system-environment interaction (i.e. Lamb shift Hamiltonian) 
have been neglected in Eq. ( |A.l| l, and that the rotating wave 
approximation is implicit its derivation. Therefore the time 
scale of the system ts ~ max {g^\ w^^' ) must be much smaller 
than the dissipation time, i.e. ts <& 7"'. 

Our dissipative system-environment coupling operators 
yfycTx^ give rise to six open decay channels (for each a), 
associated with the frequencies {+Wq,, +(jJa ± g]- Consider, 
for instance, the Lindblad operators within the cold dissi- 
pator Dc- While the operators Ac^±ai, produce transitions 
|0,,0/,0f> « |Oh.O/Jc), |1w1/,0c) « IIh-UIc), in which the 
cold bath exchanges energy locally with the cold qubit only, 
the remaining operators, e.g. Ac^±a)^.+g, are instead responsi- 
ble for processes like |1m.0/,0c) <-> (llwOftlc) ± IOu-IaOc)) / V2, 
|0,vl/,lc> (|1h'0/,1c) + |0„,1aOc)) / V2, in which bath c now 
exchanges energy with the work and hot qubit as well, in a 
delocalized way. It is in this sense, that we refer to Dci as 
modelling a delocalized dissipation effect. 

Of course, as the limit of vanishing coupling g is ap- 
proached, the rates Fq. -t„+g and FQ..-t„_g become equal, and all 
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FIG. 2: (Color online). Performance characteristic of the quantum 
absorption refrigerator, according to the localized (dashed blue) and 
delocalized (solid red) descriptions of the dissipation, characterized 
by Refs. ISHTl and Eq. l |A.l^ respectively. The operation tempera- 
tures are Ti, = 66.25 and = 4.78. We also fix T„ = 127.33, 
w,, = 56.87, g = 0.1 and y = 10"'' and p, = 10"^ The only re- 
maining free parameter, 0)^, is varied from to the upper limit of the 
cooling window (^c.max = cdAT^T,,, - Ti,)]/[T,,.(Ti, - 7^)] = oj^JcIt 
and the cooling power for each configuration is plotted versus 
the corresponding COP, normalized by the Camot bound sc- The 
cooling powers are also normalized by their maximum values Qc^^na^: 
8.39 X 10"^ and 1.27 x 10"* respectively. Clearly, uj^ = corresponds 
to <3c = £ = 0. In the ideal case of Eq. ([2j, co^ — > oJcmm results in 
£ — > Ec and — > 0. However, any irreversibility would yield a COP 
not monotonically increasing with toi^, and therefore, a multi-valued 
or even closed performance characteristic. 



delocalized transitions tend to compensate each other. For 
^ = 0, only two (localized) decay channels remain open for 
each bath, namely A(j, -t„^(0) oc cTa^^: ® !„, which stands for 
the usual ladder operators for qubit a (the remaining qubits 
are denoted as a). The idealized model of Refs. fS^-^ is then 
recovered in this limit. Note that since the frequencies appear 
exponentiated in the spectral coiTelation tensor, delocalized 
dissipation effects are intuitively expected to be still relevant 
(i.e. r-t„^+^ ^ ^±aj„-g) even for arbitrarily small g «K 1 ll25l . 

Universal performance bounds. — Equipped with the sta- 
tionary solutions of the Markovian master equation Eq. ( |A.l| l, 
we can compute the central quantities of our study, namely 
the rates at which energy from each bath is fed into the 
system, i.e., the heat cuiTents. These are given as Qa - 
Tr{7/ref2D[v ED- In particular, we refer to Qc as cool- 
ing power. Then the COP simply reads: e = QdQw HHU. 

Let us comment on our intuition linking delocalized dis- 
sipation with irreversibility in the operation of the machine. 
The Camot COP sc is realized at the upper limit of the cool- 
ing window, whenever s = Qc/Qh - <^cl<^w |6|. This holds 
in principle for arbitrary g, if Eq. is replaced with 



the qubit-qubit interaction, no matter how small, it becomes 
impossible to approach sc arbitrarily closely, as illustrated in 
Fig. |4] This suggests that the unavoidable delocalization in 
the dissipative dynamics makes the refrigerator nonideal and 
somehow wasteful, thus introducing a fundamental source of 
iiTeversibility that prevents it from cooling at the Camot COP 
in practice. 

It is therefore crucial to investigate whether, instead of the 
unreachable ec, there might exist a realistic tight bound on 
performance, such that attaining it would mark the function- 
ing of the refrigerator as effectively optimal. We believe that 
a sensible figure of merit to investigate in this context is, for 
instance, the COP e, at maximum cooling power (3c,max- One 
can then seek to devise a general upper bound for such a per- 
formance indicator, and to characterize a region within the 
space of the control parameters {w,,,, r,,,) that allows the sat- 
uration of the bound. This would provide useful workpoints 
for the efficient implementation of the machine. 

Let us first consider the ideal model of |5 -7 1 based on local- 
ized dissipation. We performed an extensive numerical anal- 
ysis and obtained a neat and tight bound for the COP at max- 
imum power, e, < ec/2 (see Fig. |5(a)| ). The bound is uni- 
versal, in the sense that all free parameters of the refriger- 
ator were picked completely at random, only respecting the 
physical constraints imposed by the Born-Markov and rotat- 
ing wave approximations. By resorting to the exact stationary 
state derived for this model in ||6l, one can prove the existence 
of the bound analytically. Namely, given r„, and w„, such that: 
(i) (Wn,/r„. /, <K 1, and (ii) <w„,/t <sc 1 (cf. caption of Fig.|4]l, 
one has that the COP at (3c,max is approximately given by 
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wich saturates to 1/2 for Tc/Th <K 1 11251 . Note that we can- 
not let Tc — > 0, as it lies outside the range of application of 
our model. Conditions (i) and (ii) can thus be viewed as de- 
sign guidelines that guarantee the best possible performance, 
in the regime of large temperature difference, and in the ideal 
case. Note as well that these are just sufficient conditions for 
optimal performance, since e.g. even the machine in Fig. [4] 
cools close to the bound although a>H./Ti, ~ 1. 

When Eq. ( |A.l| i is used instead to account consistently for 
the dissipative dynamics, which is the most interesting case, 
our extensive numerics show that e, is upper bounded by 
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g = -i [//ref, g] +Za^a® Iff [p], 



(2) which is again a universal result, independent of the details of 



where the local dissipators 



£)Sf or the ones used in 15 - 7 



la can be, for instance, our 
. Recall from the preceed- 
ing considerations, that such a localized model for the dis- 
sipation is only physically consistent in the limit of strictly 
vanishing g 11261 . On the contrary, if the realistic delocal- 
ized description of Eq. (|A.1[) is adopted, given any value of 



the system and the temperatures of the baths (see Fig. 5(b) i. 
Likewise, the bound is saturated in the large temperature dif- 
ference regime, when conditions (i) and (ii) are jointly met. 

Interestingly, despite of the iiTeversibility derived from 
delocalized dissipation, the optimal design guidelines ob- 
tained under the simpler localized description of ||5}^7| still 
remain valid. Our results also suggest that the particulars 
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FIG. 3: (Color online). COP at maximum cooling power for 10^ 
random refrigerators (a) following a localized dissipative scheme 
as in (5}Q, and (b) the delocalized one resulting from Eq. l |A.l^ . 
The operation temperatures and T;,, as well as T,,., w,,,, g, y and 
Pi = p, were chosen completely at random, always satisfying the 
constraints kgT^^ <c y"' (p^') (Bom-Markov approximation) and 
g"' <K {y P^'} (rotating wave approximation). The value of a>c 
yielding <3c,max was found in each case (therefore fixing (jJi, = cJc+cj^) 
and the corresponding e,/ec, plotted. The COP at Qc.mnji is upper 
bounded by ec/2 and 3ec/4 in these localized and delocalized dissi- 
pative models, respectively. 



of the heat baths, encoded in the spectral correlation tensor, 
may be behind the raise of the normalized e» to significantly 
higher values. The specific role of the environments in the 
enhancement of the COP at maximum power is worthy of 
deeper investigation and will be studied elsewhere. 



ing several interacting parts. Our results provide evidence that 
all these machines should inherit this fundamental limitation 
in their performance, in contrast with self-contained thermal 
devices modelled through localized dissipative schemes 1 1 -4|. 
As an alternative to sc, a more useful universal upper bound 
for e, at maximum cooling power Qc,max, was then obtained 
together with sufficient conditions to saturate it in the limit of 
large temperature diff'erence. The efficient performance of the 
studied machines was not found to relate in any obvious way 
to stationary quantum correlations present in the system. We 
emphasized how the adoption of different environmental noise 
models may raise the upper bound on the COP at maximum 
power. This intriguing phenomenon can render the practical 
realization of highly efficient quantum refrigerators more ac- 
cessible with present-day technology lfT4l[T5l and will warrant 
further investigation in the future. 
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Stationary quantum correlations. — We finally investigate 
the stationary bipartite quantum correlations established in 
the refrigerator, focusing on the realistic dissipative model 
of Eq. (A.l I. Given the structure of the reduced states 
within the 2x2 bipartitions w-h, w-c and h-c are diagonal 
and therefore, unentangled and completely classical. We then 
consider the only two-dimensional subspaces of (^^'Ks.o 
which are in direct interaction according to H^^f, namely those 
corresponding to the machine virtual qubit and the cold qubit. 
While no entanglement is found, more general quantum 
correlations measured by quantum discord ll22l l23l |271 1281 
are always present in this relevant bipartition ||251 l29l . The 
structure of stationary quantum discord, analyzed in detail in 
the Supplemental material |25 1, does not however exhibit any 
specific relationship either with the maximization of or s, 
or with the behavior of e, or (3c,max, as the control parameters 
{cjw, r„.) are varied. This supports the conclusion that the 
only essential quantum ingredient exploited by the machine 
is in fact the discreteness of its energy spectrum fSH?!- 

Conclusions. — We studied the three-qubit quantum absorp- 
tion chillers introduced in |I5}{3 adopting a physically mean- 
ingful system-bath interaction model. The resulting delocal- 
ized description of dissipation prevents the refrigerator a pri- 
ori from cooling arbitrarily closely to the Carnot COP Sc- The 
three-qubit chiller may be seen as the prototype of a whole 
family of more complex self-contained refrigerators compris- 
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SUPPLEMENTAL MATERIAL 



latter time as long as the Born approximation holds (see 
below). 

The interaction picture Liouville-Von Neumann equation 
for p it) is then suitably manipulated. Next, the Born or weak 
dissipation approximation, according to which p(t) ^ g (f) ®x 
is performed. The Markov approximation, that consists in ne- 
glecting any memory effects in the reduced evolution, finally 
leads to (see [24| for details) 



a,l3 

(A.2) 

where ^(f) = e'^'''g{t) e^'^''' stands for the interac- 
tion picture reduced state of the refrigerator qubits, and 
HD,a(t) = e'^'^'Ho^aC-'"'". The dynamical map (D(f,0) 
evolving g{0) into ^(f) that results from Eq. ( |A.2| l, 
has the semigroup property under map composition: 
O (r, 0) ■ O (i, 0) = O (f H- s, 0), which implies that Eq. \h.2) 
can be cast in the standard Lindblad form LA4J of Eq. ( |A.1[ ). 

In order to achieve this, we shall decompose the system 
operators cr^^ from the system-baths interaction term Hi - 
Y^a (^x„ ® Ba into eigen-operators J[a,ai of //ref such that 



We begin this supplemental material by systematically 
constructing the IVIarkovian master equation Eq. (2) of the 
main article. Taking as starting point the analytical formulas 
for the stationary state of the machine qubits obtained in 
Ref. in under a localized description of the dissipation, we 
also prove that the performance at maximum cooling power, 
under conditions (i) cOw/T^j, <s: 1, and (ii) liJh,/t <s; 1, is upper 
bounded by ec/2. We conclude by reporting on the observed 
properties of the stationary quantum discord between the 
machine virtual qubit and the cold qubit, as a function of 
various relevant parameters of the model. 

Derivation of the Markovian master equation Eq. (2). — 

We will now consistently build the Markovian master equation 
for the reduced state g of the three refrigerator qubits 

g = -/ [//ref,^] + 2 r„,^ iAa,cjgAl^^ - ^{Al^^Aa,cj,g}+\ ■ 

(A.l) 

The process starts by taking the interaction picture with 
respect to the free Hamiltonian Hf — ^o.« + Hj + Hsg,- 
An initial preparation, uncorrected between system and 
environment is chosen so that p (0) - g (0) ig> where 
X = ^ Xa"^ e^^"'^^" ■ This choice guarantees on the one 
hand that the reduced evolution is a completely positive (and 
trace preserving) dynamical map (CPT) [ A2| and on the other, 
that the average of the bath operator Ba = YiA sA'^a,A - <2^^) 
vanishes initially iXBBapiQ) - 0, and actually also at any 



x„ - ^-^ff.w [HieU ■!^a.oj\ - -i^^a,w (A. 3) 



The non-Hermitian Lindblad or jump operators Jia,a) ^xe, 
defined as 



where \j) is an eigenstate of Href with energy ajj and 
non-degeneracy is assumed. The eigenvalues of i/ief are 
{0, (jJu-, 2a>i,, a>c, (aJh- + oJiu ^^h + <^c, <JJh — g, <jJ/i + g], and their 
coiTesponding eigenvectors 

|1> = |0,.m), |2) = llvAO,), 13) = |l,,Ulc.>, 
|4) = |0„0,,1.>, |5> = |1„1aO,), |6) = |0,,1,,1,), 
|7) = (|l,vO,,l,)-|0„,l,,0,))/V2 
|8) = (|l,Al.> + |0,4/,0,))/V2. 

Therefore, from Eq. ( |A.4| i it is easy to see that there are 
only six open decay channels (i.e. transition frequencies oj - 
(jjj - Oik with non-zero Jia,<si) for each bath a, corresponding 
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to {+a>a, ±<j^ff ± g}- These jump operators are explicitly 

= Vr(li><2| + |6><3|) 

= Vr(|4><8|-|7)<5|)/V2 
^u..^„-, = Vr(|4><7| + |8)<5|)/V2 

^/,^„ = Vr(|2><5| + |4)<6|) 

= Vr(l7)<3| + |i><8|)/V2 
^h.c^-g^ Vr(l8><3|-|i><7|)/V2 

= Vr(li)<4| + 15)01) 

= Vr(|2><8|-|7)<6|)/V2 
= Vr(|2><7| + |8><6|)/V2. 

The remaining Lindblad operators are just given by the 
adjoint of these, since = 

The interaction picture system-environment coupling 
Hamiltonian may be now written as 

H, (f) = e'"''H,e-'"'' = ^ ® (?) , (A.5) 

where B„ (f) = Y^a gi(fl^e""^'' - fl^e'"''). Gathering all this 
back into Eq. ( |A.2[ i yields 

CC, tt>, tj' 

(A.6) 

where the (complex) bath correlations Ca.oj are defined as 

r°° 1 

Ca,oj = t/.? e'"' trB;ifB„ (f) B„ (f - .?) = -r„, <^ + /5 „ . 
Jo ^ 

(A.7) 

Note that due to our choice of the initial preparation 
and provided that the Born approximation holds, all bath 
correlations depending on Xisx^a (0 Bp (t - s) with a fi 
vanish, i.e. the baths are independent. 

The details of Ca,w may be also worked out, exploiting that 
the baths were prepared in a thermal equilibrium state ll24l . 
Choosing oc yZZ^, the spectral correlation tensor r„ reads 

r,,„cx.,V"/2(sinh^[\ (A.8) 

where the proportionality constant is of order one and may be 
absorbed into y. 

Another important step towards the derivation of the master 
equation Eq. is the assumption that the typical system 

time-scales \u) - dLt'|"' with a< (J , are much smaller than the 
relaxation time y"', which allows to discard all rapidly oscil- 
lating terms oj + u' that average to zero in a coarse-grained 



picture of the reduced dynamics (rotating-wave approxima- 
tion). This leaves (we refer again to |24| for details) 

where we have discarded the Lamb-Shift term 
~' Zo-.w'^a.w as usually done when work- 

ing in the quantum optical regime. 

The only thing that remains to be done in order to recover 
Eq. ( |A.l| i, is to transform Eq. ( |A.9| l back into de Schrodinger 
picture by noting that 

Q = -iHFe-''^'"ge''^'"+ie-'"'"ge'"'"HF+e-'^'"te'"'" . (A.IO) 

This immediatly yields Eq. ( |A.l| l if one identifies with 

p-'HfI m piHfl _ Jwt m 

The reduced dynamics generated by equation Eq. ( |A.1[ ) 
may be understood as a stochastic process in the Hilbert state 
space of the refrigerator qubits, in which the deterministic 
evolution of any pure state is interrupted by discontinuous 
quantum jumps N^^Aa,oj occurring at rates Fa oj- 

The density matrix g(t) at any time t is recovered as an 
esemble average over stochastic trajectories ll24l IA5I . 

It is clear that all interaction picture jump operators 
^a,oj introduced above, can produce delocalized dissipation 
(quantum jumps) in the sense discussed in the main article. 
As g gets closer to zero, the rates ra,±a>„+g and ra,±aj„-g lin- 
early approach each other, meaning that the jump processes 
'^a,±w+g and ^a,±aj-g bccomc equally likely, so that their 
delocalized contributions start to compensate. In the strict 
limit of g - Q, only one transition frequency (cja) is left for 
each bath, the corresponding Lindblad operator being just 
the sum Jla,±aj + ^a,±<^+g + ^a,±w-g 0^ cr„,- [cf. Eq. ( |A.4| )], 
which is a localized jump operator. Note that even if the 
difference ITaj-^^^^ - F^ j-^^^l depends linearly on g <s: 1, 
the exponentials in Eq. ( |A.8| l make it extremely sensitive 
to the qubit-qubit interaction strength, so that even a slight 
departure from the non interacting case can make delocalized 
dissipation effects very important. 

Let us finally comment on the underlying assumptions 
leading to Eq. ( |A.l| i. Even though the rotating-wave ap- 
proximation makes the problem much more tractable, it is 
not essential and one could just avoid it [cf. Eq. ( |A.6[ )]. 
Situations in which the dissipation times become comparable 
to the system time scale (i.e. the realm of quantum Brownian 
motion) may be then accounted for, conceivably resulting in 
qualitative differences. Accounting for the non negligible 
renormalization effects of the system-environment interaction 
on the system itself, becomes important in these cases. On 
the contrary, within the quantum optical regime, the rotating- 
wave approximation only slightly modifies the reduced 
dynamics but not the stationary states of the refrigerator, thus 
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leaving our results unaffected. 

If the baths are assumed to have some structure, the Markov 
assumption cannot be consistently performed. Nevertheless, 
as long as the dissipation strength remains sufficiently weak 
so that the Born approximation is still in place, no qualitative 
differences should be expected from what we report in the 
main article. Finally, if the dissipation becomes strong 
enough, the thermalization of a single isolated qubit in 
contact with its corresponding bath is no longer guaranteed, 
so that the basic operation of the refrigerator is compromised. 

Analytical derivation of the perfomance bound ec/2. — 
Considering the localized dissipative model of ||5}{7l, we 
shall now prove that its COP at maximum power e, is upper 
bounded by ec/2, whenever conditions 



TATh - T,) 

(0„<^T= — 

are met, at large temperature difference TdTi, <s; 1. 



(i) 
(ii) 



Proof. Our starting point will be Eqs. (18) and (8)-(10) in 
Ref.l6j, where the cooling power Qc was given as 



Qc,w = q- 



(A.ll) 



and 
A 

Q.ai3 



1 , g-oJalTa 1 



\+e-i"<tlTa 



1 1 



+ 



1 1 

1 1 



(A. 12a) 
(a, /Si^h) 

(ai^/3 = h) 
(A. 12b) 



Here, q, qj and Qap only depend on the three, possibly dif- 
ferent, dissipation rates while A and depend on all 
temperatures and frequencies. All we need to do is to find the 
(be that maximizes Eq. (jATTJ, and then compute the corre- 
sponding e, . 

First of all, note that conditions ^ and (|ii| imply <±>clTc ^ 1 
for any ujc < ^ax, i-e. within the cooling window 



1 > e 



-die IT , 



-«Jt.max/r,. 



^ 1 - 



(A.13) 



Due to the ordering r„, > T/, > Tc in the bath's equilibrium 
temperatures, we must also have oJc/Th ^ 1. Since w/, - 
w„ -HWc, this translates into aih/Ti, « 1. Therefore, conditions 
^ and (|ii| can be alternatively stated as 



(A. 14) 



Furthermore, in the limit of large temperature difference 
TclTh <s: 1, assuming that Tc/T/, is at least of order w„77'/, 
and a>H,/T, we also have 



CJr CL)r 
< 



(A.15) 

Th Th T Th Th 

As a consequence, Eqs. ( |A.12a| i and ( |A.12b| i can be expressed 
as 

A ^ 1 (e-- /^'. - ) + O j (A. 16a) 

a,;?- . (A. 16b) 



In this regime, the denominator of as given in Eq. ( |A.l 1[ ) 
becomes independent of aic, so that its maximization is equiv- 
alent to that of Eq. ( |A.16a[ ). This yields 



Tr 



(A. 17) 



The solution to an equation of the form jc - a may be 
expressed in terms of the Lambert-W function or product- 
logarithm IIA8L as X = W() (a e). Therefore Cjc reads 



+oj„./T„.-a>,,/T,,) 



)1 



(A.18) 



Among the properties of Wo (z), we shall make use of its series 
expansion around z = e 



Wo(z)4.^. 



(A. 19) 



Taking again into account, ^ 1 -i- lo„/T„ - 

(^k/Ti, which, combined with Eqs. ( |A.18| l and ( |A.19| l, results 
in 



(A.20) 



2 \r,, r,., 



The COP at maximum cooling power e» - Uc/Uw 1 6 1 normal- 
ized by sc, may be thus approximated by 

Ec 2ec \ Th r„, / 2 \ Thj 2 

which saturates in the limit of large temperature difference 
T.ITh « 1. ■ 

Stationary quantum correlations. — We report now on the 
quantumness of the stationary bipartite qubit correlations 
in the refrigerator, resulting from our realistic delocalized 
modelling of the dissipation. 



The dissipative dynamics of Eq. ( |A.l| l annihilates 3-qubit 
X-states that have = as the only nonzero matrix 
elements outside the main diagonal when expressed in the 
computational basis. If any of the qubits is traced out, the 
remaining 2x2 density matrix is diagonal and only involves 
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the populations of g°°. Therefore, entanglement and quantum 
discord in the bipartitions w-h, w-c and h-c vanish. 

It is interesting to look instead to the only two-dimensional 
subspaces of l^s^a which are placed in direct interaction 
through Hi, that is, the machine virtual qubit and the cold 
qubit. The corresponding reduced state, reads 

P = |1,vOa><1,AI + |0„,1a><0,,1,,| + |0,><0,| + |1,><1,| . 

When expressed in the basis {|l,t.O/,) , |0„l/,))®{|Oc) , |lc)), the 
ate is eiven bv 



wnen expressea m tne oasis i|i 
2-qubit X-state pJJ is given by 



Qvc = 



1 

TV 



^ 

eel 936 

Q'^^ £.^3 

£.^4 



(A.22) 



where N is the normalization factor According to the 
positivity-of-the-partial-transpose separability criterion 
IATOI . the state is entangled iff 



^36 > 



(A.23) 



However, our stationary states are such that <K for 
7 G {!,••• ,8}, and therefore, no bipartite qubit entanglement 
can be present in them. On the contrary, one always finds 
nonzero stationary quantum discord between the machine 
virtual qubit and the cold qubit. In the search of the least 
disturbing local measurements for its quantification, we will 
restrict to projective measurements only. Since is an 
X-state, its discord can be computed analytically, using the 
formulas of Ref. llAlli . 

In the first place, we fixed Ta and cj^^-, and looked at the 
stationary total [J(£),,f)], quantum [Digy^)} and classical 
[licTvc)} correlations between the machine virtual qubit and 
the cold qubit, to see whether they play any role in the 
maximization of the cooling power Qc for < u)c < Wc.max 
(see Fig. Similarly to Qc, for intermediate values of oj^. 
all correlations are peaked around a certain value of 0)^ that 
nevertheless usually differs from Cjc [cf. Figs. 4(c) and 4(d)| . 
Smaller work frequencies yield a monotonic behaviour of 



correlations instead, as shown in Fig. 4(b) while larger values 
of w„. reveal a more intricate structure [see Fig. |4(f)| . 



In Figs. 4(b) and 4(d)| the measurements that maximize 
I{o-vc) consist in e.g. projections onto |+c) = (|0c) + |lc)) / V2 
and l-c) = (|0c) - |lc)) / V2 for any w,.. However, in Fig.jiffl 



projective measurements in the computational basis of the 
cold qubit, become optimal in the interval 14 < < 66. 
These discontinuous changes in the optimal measurement 
schemes result in a non difFerentiable classical correlations 
and quantum discord. In all three cases, the COP increases 
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FIG. 4: (Left column) cooling power as a function of ai^ at fixed 
r„, a),,., g and y and (right column) total (dot-dashed purple), classi- 
cal (dashed red) and quantum (solid blue) correlations for T„ = 180, 
= 95, = 80, g = 0.1, 7 = 10-^ and io,„ = 10 [(a) and (b)], 
w„. = 15 [(b) and (c)] and w„, = 30 [(d) and (e)]. The gray dotted line 
marks the position of ai^. 



linearly with w^, and starts to decrease only as w^max is 
approached. It seems therefore, clear that the maximization 
of Qc and s are not related in any way to the only non 
vanishing 2x2 stationary quantum correlations in the system. 

We also ruled out any possible interplay between quantum 
discord at ajc and the maximization of e, and Qc,,mx, using 
(jj„ and Th as control parameters. In Fig. |5]we plot e, and 
the corresponding quantum discord as a function of <ij,v for 
different temperatures r„,. After increasing abruptly for small 
work frequencies, the COP at maximum power starts to decay 
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FIG. 5: (a) COP at maximum cooling power e, and (b) quantum 
discord at Hj^ as a function of a>„. for different temperatures T„.. Pa- 
rameters were chosen as T;, = 17, r^. = 13, g = 0. 1 and y = 2.5 x 10"' 
and r,, = 50 (dotted green), T,, = 100 (dot-dashed blue), r„, = 150 
(dashed orange) and r„, = 200 (solid red). 



yet build more quantumness in the system if one also leaves 
w^. as a free parameter 

From all the preceeding we see that there is no clear rela- 
tionship between the quantumness of the bipartite qubit-qubit 
correlations established in the stationary regime and the ef- 
ficient performance of the refrigerator. In the idealized case 
of vanishing g, the ability of this quantum machine to satu- 
rate the Camot COP comes from the discreteness of its energy 
spectrum, in contrast with any continuous "classical" counter- 
part Q. Actually this discreteness and thermalization are the 
only essential building blocks for the basic operation of the 
refrigerator It may not be surprising then, that any residual 
quantumness of correlations appears in the system as a by- 
product rather than as fundamental resource for its enhanced 
performance. 



as oj„ grows. When it comes to its temperature dependence, 
e* seems to increase with Th, as shown in Fig. |5(a)| In 
Fig. 5(b) we see that for small Wh,, quantum discord is also 
an increasing function of the work frequency. The optimal 
measurement scheme changes from {|+c) (+cL (-d) 
to {|0c)(0c|, |le)(lfl) at some uj^, which produces a sharp 
maximum. For > w,,,, discord decays monotonically. At 
any fixed oj^y, it is decreasing with the work temperature. The 
corresponding <3c,,nax grows exponentially with a>„ and also 
increases as a function of r„. 

Even if the maximum discord at fixed r„, does not coincide 
with the maximum of e», it still marks a useful operation 
point of the refrigerator where a certain compromise between 
<3e,max ^nd is achieved. Our results also suggest that the 
COP at maximum power and the corresponding cooling 
power are enhanced, at fixed 0)^. (and sufficiently small g), at 
the expense of the destruction of quantum coiTelations. It is 
possible, however, to increase these two figures of merit and 
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